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I.  INTRODUCTION 

Our  previous  work,  hereafter  referred  to  as  Report  II,  concerned  the  beam 
attenuation  due  to  "seif-stripping"  resulting  from  an  "optimum"  electric  field 
having  the  components  due  to  transverse  and  longitudinal  polarization, 

Ey  -  vb  B  sin6  (D 

E2  *  vb  B  cos28  sin6  (2) 

Ex  *  -vb  B  cos  6  sin2e  (3) 

where  vb  is  the  beam  velocity,  and  B  ■  /bx2+b22  ,  b2  ■  B  sin6 ,  Bx  *  B  cos6 
(See  Figure  1). 


We  assumed  a  confinement 
time  T|  of  the  form 
2Rn 


T| 


vb  sin^e 
with  e  <  90*,  and  R^,  being 
the  beam  radius. 


(4) 


Figure  1 


We  also  assumed  another  confinement  time  T  given  by 

£ 

cos (wt )  *  1  - ,  for  6  <  90°  (5) 

K 


gB 

where  w  ■  —  and  R  are  respectively  the  electron  gyrof requency  and  gyroradius. 
“e 

In  the  present  work,  we  will  study  alternative  models  for  the  electric 
field,  and  the  corresponding  confinement  times.  Our  purpose  is  twofold. 

First,  we  wish  to  study  the  effect  of  different  reasonable  models  on  the 
confinement  time.  Second,  we  wish  to  combine  the  confinement  time  (T)  with 
the  (energy  dependent)  collision  cross  section  (o)  to  obtain  an  estimate  for 


1 


T  T 

the  total  beam-induced  stripping  probability  P  * -  -  — —  as  a  function  of 

onv  oj 

the  beam  parameters. 

II.  ELECTRIC  FIELD  KITH  TRANSVERSE  COMPONENT  ONLY:  CONFINEMENT  TIME 

In  the  first  model,  the  electric  field  has  only  one  component,  due  to 
transverse  polarization.  More  specifically,  we  have 
Ey  -  fvb  B  sin6,  CKf<l 
Ez  -  0 

Ex  -  0  (6) 

where  f  can  be  regarded  as  field  strength  parameter,  f  *  1  corresponds  to  a 
fully  developed  polarization  field,  while  f<l  describes  reduction  due  either 
to  ambient  plasma  screening  or  to  transient  situations  where  the  electric 
field  has  not  fully  developed. 

The  resulting  electron  velocity  equation  becomes 


x(t)  *  (1-f)  vb  sin^Q  •co8(uit)  +  vb  (fsinS'cosS)  (7) 

y(t)  •  (1-f)  vb  sin0  •  sin(u)t )  (8) 

z(t)  ■  (1-f)  vb  sin0*cos0  Ll-cos(u>t)j  (9) 

For  the  field-free  case,  f  -  0,  we  have 

xQ(t)  ■  vb  sin^0  *cos(uit)  ♦  vb  cos^0  (10) 

y„(t)  ■  vb  sin0*8in(wt)  (11) 

z0(t)  *  vb  sin6*cos0  [l-cos(u>t)j  (12) 


•  • 

For  the  "maximum"  strength  field:  f  ■  1,  we  have  x(t)  ■  vb,  y(t)  *  0,  and 
z(t)  •  0  which  describes  undisturbed  electron  motion. 


Integrations  of  Eqs.  (8)  and  (9)  give  the  trajectory  equations 

y(t)  “  (1-f)  R  sinfi  [ l-cos(wt)]  (13) 

z(t)  ■  (1-f)  R  sin8«cos8  [cue  -  sin(u)t)],  (14) 

where  again  u>  and  R  are  respectively  the  electron  gyrof requency  and 
gyroradius. 

The  y(t)  is  sinusoidal  and  confined  between  0  and  2(1 — f )  R  sin6,  while 
z(t)  has  a  secular  behavior  due  to  drift. 

The  confinement  time  can  be  considered  in  the  following  manner.  Looking 
along  the  positive  x-axis,  we  can  interpret  the  interception  of  y-z  trajectory 
with  the  beam  (radius  R0)  circumference  as  equivalent  to  confinement  time  Tf 
(See  Figure  2)  u 


[  “fej/ 


Figure  2 


Ue  can  now  write 
Y2(Tf)  +  Z2(Tf)  -  R02 
or  more  specifically 

2 

sin20  [  l-co8(u)Tf )]  2  +  sin20 ’cos20  [ujt  -  sin(u)Tf)]2  ■  (i  £)2  ’  ^ 

where  Y  ■  R0/R.  Computer-generated  results  from  Eq.  (15)  in  the  form  of 
Tf  *  Tf (0 ,Y )  are  given  in  Figs.  3(a)  to  3(n).  Assuming  Rc  to  have  a  fixed 

*  2 

value  of  10cm,  we  have  Y  ■  — - - ,  with  n  being  the  beam  energy  E^ 

in  unit  of  MeV. 


irw-ir*  hr+w".w*.v*  ir*ir*  »r»  «n*mm 


It  is  worth  noting  that  T0  represents  the  confinement  time  for  the  zero- 
field  case,  while  Tf  represents  the  confinement  time  for  “optimum” 
y-component~only  field;  Tf  is  infinite ,  since  the  optimum  y-component-only 
field  results  in  undisturbed  electron  drift. 

However,  Tf  can  be  considered  in  another  way.  From  Eqs.  (13)  and  (1A), 
we  write 

Y  •  (1-f)  R  sin6  (1-cos  wTf )  ■  |  (l-coswTf), 

Z  ■  (l~f)  R  sin6*cos0  (a)Tf-sinu)Tf )  ■  (wTf  -  sinwTf), 

where  ^  «  (1-f)  R  sin0  and  ^  -  (1— f )  R  sin0*co60.  Rearranging  terms 
and  squaring,  we  obtain  a  quadratic  equation  for  Tf, 
w2Tf2  -  2wBZTf  +  Tf  +  B2Z2  -  2AY  -  0, 
which  yields  the  solution 

Tf  -  ^  +  [2AY  -  A2Y2)1/2  (16) 

It  i6  instructive  to  compare  Eq.  (16)  with  Eq.  (4).  If  we  let  Y  =  0,  we  get 

,  BZ  BRq  2Rq  T| 

Tf(Y-O)  -  —  -  —  “  d_f )VbSln2e  (1-f )»  (17) 

where  we  have  made  use  of  the  relation, 

Z2  +  Y2  -  Ro2,  i.e.,  Z  -  R0  for  Y  =  0. 

For  the  zero-field  case(f  ■  0)  we  obtain 

Ty  -  To(Y«0) 

Thus,  we  have  clarified  the  approximation  Involved  in  the  confinement  time  Tt 
as  used  in  Report  II.  In  other  words,  Ty  is  equal  to  the  zero-field 
confinement  time  with  the  approximation  of  having  Y  ■  0.  In  the  following  we 
will  examine  closely  various  types  of  confinement  times  and  their  comparative 


merits. 


The  switching-on  of  optimum  field  has  the  effect  of  rotat ing  the  y-z 
trajectory  about  the  x-axis  through  90°  and  reducing  z  projection  by  a  factor 
of  cos26.  Specifically,  we  have 

r02(t)  -  y02(t)  +  z02(t)  (23) 

re2(t)  -  ye2(t)  +  zG2(t) 

■  cos^8  y02(t)  +  z02(t)  (24) 

which  shows  that  re(t)  is  less  than  r0(t)  except  at  6  -  0°.  Consequently,  the 
confinement  time  Te  from  re(Te)=R  is  longer  than  the  field-free  value  T0  from 

r0(T0) “R0 * 

To  obtain  Te ,  we  set 

Ye2(Te)  -  ye2(Te )  +  zG2(Te)  -  R02  (25) 

which  can  be  re-expressed  as 
v2 

— ^ - r—  *  u)2Tf2  -  2ujT£ ‘sinwlc- •  -  2cos8 * coswTe 

sin28  coz26 

+  2co628  +  sin28  •  sin2inTL  (26) 

R0  ^ 

where  the  dimensionless  parameter  y  *  —  *  6.34517  x  10  »n,  with  Rq  *  10cm 

and  n  is  beam  energy  in  unit  of  Mev. 

Computer-generated  results  for  Te  are  shown  in  Figs.  5(a)  through  5(h). 

It  is  worth  noting  that  at  8  *  45®,  Te  is  about  1.2xl0-7  sec  for  beam  energies 
ranging  from  1  to  70  Mev. 

IV.  ELECTRIC  FIELD  WITH  FRACTIONAL  STRENGTH:  IMPACT  ENERGY. 

A  field  of  less-than-optimum  strength  can  be  expressed  as 

Ey  “  f  vjj  B  sin8  (27) 


Ez  -  gvb  B  sin8,cos28 


(28) 


(29) 


Ex  -  -gvb  B  sln20*cos0 

where  both  f  and  g  are  respectively  0<f <1  and  0<g<l.  f  and  g  can  be 
interpreted  as  before  in  Section  II,  as  representing  ambient  plasma  screening 
and  transient  effects. 

The  resulting  (electron)  velocity  equations  are 

x(t)  •  -dvbsin20  •[  l-cos(u>t)  ]  +  gvb  sin20  •  cos@  •  sin(u>t )  +  vb  (30) 

y(t)  *  6vb  sin© •sin(ott)  +  gvb  sin0*cos0  [l-cos(wt)]  (31) 

•  _ 

z(t)  ■  6vbsin0*cos6  [l-cos(a)t)]  -  gvb  sin9*cos29  •sin(wt)  (32) 

where  we  have  written  6  *  1-f, 

Integration  of  Eqs.  (31)  and  (32)  gives 

y(t)  ■  6R  sin0  •[  1-cos  (wt )  ]  +  gR  sin0  *coz9  •  [tut-sin(u>t )  1  (33) 

z(t)  ■  6R  sin0  *cos0  •  Iwt-sin(wt)  ]  -  gR  sin0  *cos29  [ l-cos(u)t )  ]  (34) 

vb 

where  we  have  again  written  —  ■  R,  the  electron  gyroradius. 

Now  both  y  and  z  have  secular  drift  components.  The  relative  electron- 
beam  velocity  components  are 

vx  ■  -6 vbsin29 [ l-cos(wt ) )  +  gvb  sin20,cos0,sin(u>t)  (35) 

Vy  ■  6vbsin6  *sin(u)t)  +  gvb  sin9  *cos0  •[  1-cos  (u>t )  ]  (36) 

vz  “  6 vbsin6 .cos6 •  (l-cos(a)t)]  -  gvb  sin9'co620 •sin(wt)  (37) 

Taking  the  time-average  of  the  above  quantities  over  the  confinement  time 
T,  we  have  the  time-averaged  relative  electron-beam  velocity  components: 

<vx>  “  ~dvbsin20.[l  -  ~  sin(u)T)]  +  gvb  sin20  ~  (l-cos(wT)]  (38) 

<Vy>  ■  6vbsin29*  (l-cos(wT)]  +  gvbsin0  *cos0  •  [  1  -  ~  sin(u>T)l  (39) 
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(40) 

<Vjj>  *  6vbsln6 ‘cose  •  (1  -  ~  sin(iuT)]  -  gvbsin6 • cos20 •  —  [1-co6(wT)] 

As  in  our  preceeding  work,  Report  II,  we  define  the  electron-beam  "impact 
energy"  Ej  as 

Ej  ■  -|«<v>2  •  -|id  [ <vx>2  +  <vy>2  +  <Vj>2]  (41) 

From  Eqs.  (36),  (39)  and  (40),  we  obtain 

EI  "  Eb[  (1-f )2sin26  +  g2sin26  *cos26]f(T)  (42) 

1  -> 

where  Eb  ■  y  mvb4  is  the  beam  energy,  1-f  -  6  and 

F(T)  -  2(^)2[  l-cos(o)T)j  -  2(~)sin(wT)  +  1  (43) 

More  specifically,  we  rewrite  Eq.  (42)  for  three  types  of  electric  field 
models,  namely,  zero-field,  y-component-only  field  and  optimum  field. 

(i)  Zero-Field: 

We  have  f  ■  0  and  g  ■  0 

EI0  *  Eb  s1"26  F(T0)  (44) 

(ii)  Optimum  Field: 

We  have  f  ■  1  and  g  *  1 

EU  •  Eb  sin28  •cos20  *F(TC  )  (43) 

(ili)  Y-Component-Only  Field: 

We  have  g  ■  0, 

Ely  -  Eb(l-f )2-sin20*F(Tf)  (46) 
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These  three  equations  represent  the  field  models  with  their  corresponding 
confinement  times,  and  will  be  considered  in  detail  in  the  following  section. 

V.  COMPARISON  OF  FIELD  MODELS  AND  CONFINEMENT  TIMES 

In  Report  II,  the  zero-field  model  is  represented  by  Eq.  (40)  of  Report 
II,  i.e., 

<v  >2 

— — — —  ■  sin20{ 2(^—  )2[  l-cos(^e)]  +  1—2 —  )sin(i|>e)} 

vb2  Ve 

which  can  be  rewritten  as 
EI0  -  Ebsin20  * F ( Tu ) 

1  o  1  O 

where  Ejq  “  ^“e^eb^  o»  "  ]  mevb  •  anc*  ^e  *  wTn  •  We  have  now  used  the 

correct  confinement  T0  (as  given  by  Fig.  6)  instead  of  the  approximate 

2R0 

confinement  time  Tj  ■  • 

Let  us  now  compare  the  two  confinement  times  T#  and  T0.  As  pointed  out 
in  Report  II,  T|(  is  symmetric  with  respect  to  0  *  45°,  but  is  not  valid  for 
6  =  90°.  Both  Tu  and  T0  are  infinite  as  0  •  0°  as  they  should,  however,  T0 
approaches  asymptotically  its  lowest  value  at  0  ■  90°.  Note  that  T0  is 
obtained  from  Eq.  (15)  with  f  ■  0,  and  re-entry  of  electron  into  the  beam  is 
not  considered.  Therefore,  at  0  =  90°  and  especially  for  R  <  R0,  T0  is,  like 
T ,  not  valid. 

For  numerical  comparison,  let  R0  ■  10cm,  Eb  ■  1  MeV  and  0  *  45°;  we  then 
get  T0  *  4.45xl0“®6  and  T|  «  1.44xl0“®s.  Note  that  T0  >  T||  for  0  <  80°. 

The  result  for  the  optimum  field  model  in  Report  II  was  given  by  Eq.  (44) 
of  Report  II,  which  is  rewitten  in  the  form 
Ele  -  Ebsin20*cos28  F(T[ ). 


.V.V. 
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Again,  in  our  present  result,  as  given  by  Eq.(44), 

Ele  *  Ebsin28 'cos^  F^Te^« 

we  have  used  the  correct  confinement  time  Te  instead  of  the  approximate  T|. 

Figs.  5(a)  through  5(h)  show  that  Te  is  almost  constant  with  a  value  of 
about  1.25xlO~7s  for  8  ranging  from  45°  to  75®  and  beam  energy  ranging  from 
1  MeV  to  70  MeV.  Te  increases  rapidly  as  0  approaches  0°  and  90°.  For 
numerical  comparison,  we  let  Rq  -  10cm,  Eb  -  1  MeV  and  6  «  45°;  we  then  get 
Te  *  1.25xlO“7s,  thus  we  have  Te  >  T0  >  Tj . 

Now  we  come  to  the  interesting  and  importand  result  as  represented  by  Eq. 

(45), 

Ely  -  Eb[ (l-f)2sln20j  F(Tf ). 

Note  that  the  "strength"  parameter  f  plays  an  important  role  in  determining 
the  value  of  Eiy.  Especially  for  f  >  0.9,  Tf  increases  drastically  and,  as 
expected,  la  infinite  at  f  -  1.  [See  Figs.  3(a)  through  3(n)]. 

Again  for  numerical  comparison,  we  let  R0  •  10cm,  Eb  *  1  MeV  and  8  -  45°; 
we  then  get  Tf-o.9  =  1*6  x  10_7s,  thus  we  have  Tf.Q.9  >  Tt  >  To  >  T |, . 

Comparing  the  impact  energies  as  given  by  Eqs.(44),  (45)  and  (46),  we 
have,  for  R0  »  10cm,  Eb  ■  1  MeV  and  8  *  45® 

EIo  -  E5  sin26*F(T0)  *  (1  MeV)(7^=)2(0.2)  *  100  KeV, 

Ei£  ■  Eb*sin26 *cos28 >F(TC ) 

«  (1  MeV)«(/y-)2-(0.48)  *  120  KeV 
Eiy  -  Eb(l-f)2-sin28-F(Tf) 

*  (1  MeV)*  ( l-0.9)2(^»-)2(0. 3)  *  1.5  KeV 
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It  is  worth  noting  that  while  the  relative  impact  energies  Eiq  (for 
zero-field)  and  Eje  (for  optimum  field)  are  comparable,  Ejy  (for 
y-component-only  field  with  f  ■  0.9)  is  two  orders  of  magnitudes  smaller. 
This  significant  quantitative  difference  between  Ejy  and  both  Ejq  and  Eje  is 
brought  about  by  the  strength  parameter  f,  and  we  will  discuss  the  various 
consequences  resulting  from  the  differences  among  the  three  impact  energies 
Ely*  EI0  ant*  Ele* 

VI.  DISCUSSION  AND  CONCLUSION 

We  start  with  the  confinement-time-dependent  function  F(T)  as  defined  by 
Eq.  (43), 

F(T)  -  2(^)2  [  l-cos(u)T)]  -2(~)sin(«T)  +  1, 

where  the  confinement  time  T  is  set  equal  to  T0,  Tf  or  TE'in  accordance  with 
the  respective  field  model  [in  Report  II,  T  was  set  equal  to  Tu ] . 


Writing  ^  ■  u>T,  we  get 

F(*)  -  2(^-)2  ( 1 — cos^ )  -2(-^)siniJ<  +  I, 


(43') 


which  i6  graphically  represented  by  Figs.  7(a)  and  (b).  Note  that  F(^)  is  an 
oscillating  function  with  decreasing  amplitude  and  approaches  unity  as  t  +  "• 
It  has  a  maximum  value  of  1.6  at  ip  ■  4.2,  and  approaches  zero  as  ■  0. 

In  terms  of  the  confinement  time  T,  we  have  F(T)  ♦  1  as  T  ♦  •.  Since  we 
have  assumed  a  value  of  5x10“^  Te6la  for  the  geomagnetic  field  throughout  our 
work,  we  get  w  -  8.78x10^  rad/sec.  Thus,  F  is  maximum  at  T  ■  4.78xl0-7  sec. 
For  ♦  *  0.0878  or  T  *  10“®  sec,  we  have  the  following  approximate  version  of 


Eq.  (43* ): 


F(T )  *  -~-(ta)T)2 


(43") 


BHWHWUMMWM!*lwfcl,UH”1""** 


Rewriting  Eqs.  (44)  and  (45)  of  Report  II,  together  with  Eqs.  (44),  (45) 


(46),  we  have 

E*10  “  Eb  ain28‘F(T|  ) 

(44*) 

E10  -  Eb  ain28*F(T0) 

(44) 

Ei  *  Eb  sin20*cos20  F(T|j  ) 

(45’) 

Eie  »  Eb  sin^B'cos^e  F(Te) 

(45) 

Eiy  -  Eb  (1-f )2*sin26*  F(Tf) 

(46) 

For  the  zero-field  and  the  optimum  field  models,  we  have  used  the 
corresponding  "correct"  confinement  times  T0  and  TE  instead  of the  approximate 
confinement  time  T|  as  was  done  in  Report  II. 

As  shown  by  Fig.  (6),  T0  approaches  a  lowest  value  at  0  *  90°,  and 
increases  at  decreasing  8  and  Is  infinite  at  8  ■  0°  while  Tb  has  a  minimum  at 
0  ■  45°,  is  symmetric  with  respect  to  B  «  45°,  and  diverges  at  8  «  0°  (90°). 
Figs  5(a)  to  5(h)  show  that  Te  diverges  at  both  ends,  is  approximately  linear 
with  slight ly  lower  value  at  lower  range  of  8,  and  can  be  regarded  as 
qualitatively  similar  to  Tk  . 

The  Interesting  new  result  obtains  from  the  y-component-only  field  model 
as  represented  by  Eq.  (46)  and  is  discussed  below. 

First  of  all,  it  should  be  noted  that  all  the  confinement  times,  T ( ,  T0, 
Te  and  Tf  are  not  rigorously  valid  at  8  *  90°,  because  re-entry  of  electron 
into  the  beam  path  is  not  taken  into  account.  Therefore,  if  we  exclude  large 
6,  all  the  T'b  are  qualitatively  similar,  even  though  we  have  Tf  (with 
f  >  0.9)  >  Te  >  T0  >  Tk .  It  is  worth  mentioning  that  all  the  T's  decrease 
very  slowly  with  energy  above  5  Mev. 

Because  of  factor  ( 1 — £ ) ^  in  Eq.  (46),  Ejy  is  at  least  two  orders  of 
magnitude  smaller  than  the  other  confinement  times  if  f  is  assumed  to  be  >  0.9. 
Therefore,  Ejy  results  in  a  comparatively  larger  ionization  cross  section,  as 
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Indicated  by  the  energy  versu6  cross  section  curve  [ See  Fig.  8  of  Report  II.] 
Especially  at  small  0  (6  <  10°),  we  have  larger  Tf ,  and  also  larger  o  because 
of  factor  sin^fi  in  Ejy. 

Regarding  Tf  and  o  (respectively,  the  collision  time  and  collision  cross 
section),  we  can  offer  a  semi-quantative  argument  in  favor  of  the 
y-component-only  field  model  in  terms  of  the  relation 

T  «-l - A- 

onv  o  J 

where  t  is  the  time  between  collisions,  n  is  concentration  per  cm\  v  is  the 
beam  velocity,  and  nv  •  J  is  the  beam  current  density. 

With  J  =  lA/cm^,  6  =>  10~^cm^,  we  have  t  *  10~^sec,  which  is  a  rather 
long  time.  However,  for  f  >  0.95,  6  <  10®,  and  especially  with  Ej,  <  1  MeV, 
we  obtain  a  confinement  time  Tf  comparable  to  T. 

Another  physical  consideration,  also  in  favor  of  the  y-component-only 
field  model  is  the  following.  The  time  6cale  of  beam  operation  is  probably  too 
6hort  for  the  build-up  of  "optimum"  space-change  field  as  described  in  Report 
II;  thus,  a  partial  (f  <  1)  y-component-only  field  is  a  more  likely 
possibility.  In  fact,  y-component-only  polarization  field  was  suggested  and 
observed  [see  Ref.  3  of  Report  II. j 

In  conclusion,  we  can  say  that  our  present  work  suggests  that  the  partial 
y-component-only  field  model  is  a  physically  more  reasonable,  and  thus,  is  a 
more  likely  physical  mechanism  involved  in  the  6elf-stripping.  Conversely,  if 
the  actual  polarization  field  corresponds  to  the  partial  y-component-only 
field  model,  as  studied  here,  we  can  expect  the  self-stripping  process  to  be 
significantly  enhanced. 


FIGURE  3(a) 
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